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Abstract 

Exact solutions of a family of Heisenberg-Ising spin-lattice models for a 
coupled barotropic flow - massive rotating sphere system under microcanon- 
ical constraint on relative enstrophy is obtained by the method of spherical 
constraint. Phase transitions representative of Bose-Einstein condensation 
in which highly ordered super and sub-rotating states self-organize from ran- 
dom initial vorticity states are calculated exactly and related to three key 
parameters - spin of sphere, kinetic energy of the barotropic flow which is 
specifled by the inverse temperature and amount of relative enstrophy which 
is held fixed. Angular momentum of the barotropic fluid relative to the ro- 
tating frame of the infinitely massive sphere is the main order parameter in 
this statistical mechanics problem — it is not constrained either canonically 
nor microcanonically as coupling between the fiuid and the rotating sphere 
by a complex torque is responsible for its change. This coupling and ex- 
change of angular momentum is a necessary condition for condensation in 
this spin-lattice system. There is no low temperature defects in this model - 
the partition function is calculated in closed form for all positive and negative 
temperatures. Also note-worthy is the fact that this statistical equilibrium 
model is not a mean field model and can be extended to treat fiuctuations if 
required in more complex coupled fiows. 

1 Introduction 



Consider the system consisting of a rotating high density rigid sphere of 
radius R, enveloped by a thin shell of barotropic (non-divergent) fiuid. The 
barotropic fiow is assumed to be inviscid, apart from an ability to exchange 
angular momentum and energy with the heavy solid sphere. In addition 
we assume that the fiuid is in radiation balance and there is no net energy 
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gain or loss from insolation. This provides a crude model of the complex 
planet - atmosphere interactions, including the enigmatic torque mechanism 
responsible for the phenomenon of atmospheric super-rotation. 

We will build an equilibrium statistical mechanics model to study the 
relaxation of this complex phenomenon, possibly in terms of phase transitions 
that are dependent on a few key parameters in the problem. For a problem 
concerning super-rotation on a spherical surface there is little doubt that one 
of the key parameters is angular momentum of the fluid. The total angular 
momentum of the fluid and solid sphere is a conserved quantity. We could 
either model the problem as such, that is, the angular momentum of the 
sphere is large but finite and thus can vary, or consider the sphere to have 
infinite angular momentum, in which case, it serves as a reservoir of angular 
momentum and the active part of the model is just the fluid. 

It is clear that a quasi - 2d geophysical relaxation problem will involve 
energy and enstrophy. The total energy of the fluid and sphere is conserved 
in any frame, both rotating and inertial ones. It consists entirely of the 
kinetic energy of barotropic flow plus that of the solid sphere because we have 
assumed the sphere to be a rigid solid that does not deform, and there is no 
gravitational potential energy in the fluid since it has uniform thickness and 
density, and its upper surface is a rigid lid. Conservation of relative enstrophy 
is treated here as a microcanonical constraint, modifying the classical energy- 
enstrophy theories [H] in substantial ways, chief amongst them being removal 
of the Gaussian low temperature defect while retaining the exact solvability 
of the model. 

Higher vorticity moments are considered to be less significant than en- 
strophy in statistical equilibrium models of quasi-2d geophysical flows jH]. 
A detailed variational analysis of this topic is available in [8]. 

The results in this paper is motivated by the variational analysis reported 
in [12j and the careful and detailed simulation results first reported in the 
paper [7] and later in the book by Lim and Nebus, [6]. They agree in large 
part also with the mean field theories reported in [3] and [1]. 
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2 Energy and angular momentum in the fluid- 
sphere system 

The time dependent kinetic energy in the inertial frame {Q = 0) of the 
barotropic fluid component of the above two components system is given by 



where u and v are the zonal and meridional components of the fluid velocity 
in the inertial frame. 

This kinetic energy of the fluid component in the inertial frame can also 
be written in terms of an arbitrary rotating frame as 



where Ur — u — Up is the relative zonal velocity, ip is the stream function for 
the relative flow and 



is the vorticity in the rest frame in terms of the relative vorticity lo and the 
planetary vorticity 21] cos ^; here 9 denotes co-latitude on the unit sphere S^. 
Clearly, the value of 



docs not depend on the choice of fi, except that, as will be shown next, by 
choosing Q > 0, we can conveniently measure the varying amount of angu- 
lar momentum in the fluid. However, unlike the situation for the standard 
BVE (which will be discussed next), there is clearly a difference between the 
> and Q = expressions for the rest frame kinetic energy of the same 
generalized barotropic flow. 




h[q] 




q = 00 + 2^1 cos 6 
= Ai/j + 2n cos 9 



H[q] = Ho[q] 
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Dropping the last term which is a constant and rewriting the rest frame 
kinetic energy of the fluid (in terms of relative zonal velocity u^- and merid- 
ional V — Vr), 



we observe that the second term is the projection of relative velocity onto the 
velocity of a spherical shell rotating at angular velocity Q, which is propor- 
tional to the net angular momentum of the relative flow. This term cannot be 
zero for all time since there is no distinguished value for the spin rate of the 
solid sphere when angular momentum is exchanged between the fluid and the 
sphere. In other words, after choosing some convenient fixed spin rate Q > 
to index a rotating frame, the fluid could gain or loose angular momentum 
to the sphere and this shows up in the time- varying inner product 



To continue with the calculation of energy and momentum in the two 
components system, we note that the rest frame kinetic energy of the rigid 
sphere is easily calculated to be 



where A is a constant that depends on the radius R and the density of the 
sphere and Qs is its changeable angular velocity. Angular momentum of the 
sphere in the rest frame is directly related to its kinetic energy and given by 
the linear expression 



By conservation of total kinetic energy (sum of fluid and solid sphere 
energy) in the rest frame, the sum 



even as both terms in the sum changes as Q.s changes. Similarly by conser- 
vation of total angular momentum, the sum 





An 



'S 



H[q\ + AD,s — const 




dx UrUp + BQg — const 



even as both its terms may change over time with fig. 
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2.1 The standard BVE 



For pedagogical purposes, we now compare the kinetic energy and angular 
momentum expressions for the standard barotropic vorticity model in which 

there is neither exchange of energy nor momentum with the sphere (the fluid 
component is energetically and torque-wise isolated). This is not the subject 
of this paper. The kinetic energy of the fluid in the rest (inertial) frame is 
given by 



Ho = - dx{u' + v') 

which is superficially the same expression as Hq above but now H'q is a 
constant in time. In a frame that is rotating at angular velocity the 
kinetic energy of the fluid is given by 

H'[q] = ^ I dx {u^^-\- vl) + I dx UrUp + - / dx ui 
2 752 752 2 J52 

which is again superficially similar to the expression H[q] but differs from 

that earlier expression because it is now a constant. Unlike the previous 
situation, the net angular momentum of the fiuid relative to this rotating 
frame is fixed in time and given by 



Js 



= dx UrUp — const. 
52 



Moreover, there is now a special choice of frame angular velocity (or 
gauge) for which 



= dx UrUp — 

for all time, namely the angular velocity fl' of a rigidly rotating spherical 
shell whose angular momentum equals that of the fluid. With this choice Q', 
the expression H'[q\ becomes 

= o / dx {ul + vl) + \ I dx ul 
2 J52 2 Js2 

which apart from the constant term | J^2 dx u^ has the same form as the 
rest frame kinetic energy H^. This is the often stated line that the kinetic 
energy expression of the standard barotropic vorticity model (BVE which 
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has constant energy and momentum) has the same form in any frame, both 
inertial and rotating ones. Although the kinetic energies H'[q\ = ^Jg2dx 
{u^ + v1) and H'q = ^ J^j dx (u^ + f ^) of the same fluid flow, have the same 
form, the relative velocity {ur,Vr) in the special frame labeled by Q' has 
zero net angular momentum Mq = but the same flow's absolute velocity 
{u, v) in the inertial frame has fixed nonzero angular momentum. This fact is 
properly reflected in the statistical equilibrium models for the standard BVE 
by Frederiksen et al [2] where the conservation of fluid angular momentum 
is imposed as an additional microcanonical constraint Mq = 0. 

3 Statistical mechanics and enstrophy 

To construct a statistical equilibrium model for this first system, we should 
use a formulation that is microcanonical in both the total rest frame kinetic 
energy and the total rest frame angular momentum. This formulation allows 
kinetic energy and angular momentum to be exchanged between the two 
subsystems in the relaxation process towards statistical equilibrium. How- 
ever, such a doubly microcanonical statistical ensemble is very cumbersome 
to solve. 

We therefore assume that the sphere has infinite mass and a fixed angu- 
lar velocity Q, and thus, acts as two related infinite reservoirs of rest frame 
kinetic energy H[q\ and angular momentum Mq for the fluid. This simple 
step is justified in the study of most planetary atmospheres by the relatively 
massive planetary spheres in the problem. It leads to a significant reduction 
of technical difficulties because the equilibrium statistical mechanics is now 
based on a doubly canonical ensemble in kinetic energy H[q] and angular mo- 
mentum Mq with corresponding Lagrange Multipliers or chemical potentials 
P and a. 

We observe that the key expression H[q] in this formulation for general- 
ized barotropic flows (that exchanges energy and momentum with an infinite 
reservoir) is independent of the choice of f2 > 0, precisely because it is the 
rest frame kinetic energy of the fluid. Thus, we should choose f2 > to be the 
fixed angular velocity of the massive solid sphere. For this choice of Q > 0, 
the fluctuations of Mq measure the amount of super (resp sub-) rotation in 
the fluid relative to the frame in which the solid sphere is fixed. Unlike the 
standard BVE, there is really no special choice fl' for which the net angular 
momentum term Mq/ vanish for all time. 
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This problem is still not well-posed because without fixing or limiting the 
size of the relative fiow in some suitable norm (such as relative enstrophy 
in the frame rotating), the energy H[q\ and angular momentum can in 
principle become unbounded at fixed reservoir temperatures. So we impose 
the condition of fixed relative enstrophy since it is the square of the L2{S'^) 
norm of the relative vorticity u. The classical models based on subjecting 
enstrophy to a canonical constraint leads to Gaussian models which are not 
defined at small absolute values of the temperature pT] . 

We note that fixing the relative enstrophy by a microcanonical constraint 
means two things: (1) it does not mean that the energy H[q\ and angular 
momentum Mq are fixed, and (2) it does not mean that the resulting mixed 
ensemble is intractable although in general microcanonical constraints give 
rise to great analytical difficulties. With one more minor adjustment, we 
show in this paper that the resulting models are exactly solvable spherical 
Ising and Heisenberg models that can be solved in closed form by the method 
of steepest descent. 

The last adjustment we make is to couple the two remaining reservoirs 
into one, namely fix the statistical temperature T = of a single energy 
reservoir instead of having two separate inverse temperatures /3 and a for 
the energy and angular momentum respectively. This simple step is justified 
in the following argument. Expression H[q] for the rest frame kinetic energy 
shows that net angular momentum is essentially the second and independent 
part of H[q]; the first part of H[q] is the relative kinetic energy in the rotating 
frame. That is, even as the whole H[q] fiuctuates in relaxation with respect 
to the infinite energy reservoir, the two active parts of H[q], namely the rel- 
ative kinetic energy term (1st term) and the angular momentum part (2nd 
term) exchange energy constantly even after equilibrium is reached. Thus, 
the collapse of two reservoirs into a single energy reservoir in this particu- 
lar problem, retains the physically important and statistically independent 
mechanism of angular momentum fiuctuations. 

To summarize this elementary but important material, we note that the 
kinetic energy expression used in the derivation of the spin-lattice models 
in this paper is just the changeable rest frame kinetic energy of the fiuid 
(written with respect to a frame rotating at fixed > 0) minus the constant 
term | dx Up, 



H[q] 
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In general both terms fluctuate independently as the sum H[q\ changes in 
time due to energy and angular momentum exchanges with the coupled in- 
finite reservoir. By mapping spins to local vorticity, the first term | J^a dx 
{u^ + vl) by itself gives rise to long range Ising type spin-lattice models 
without an external field. Using moreover, the analogy between magnetic 
moments and angular momentum, it is easy to see that the second term 
= dx u.fUp which represents the changeable net angular momentum 
of the fluid (relative to the frame rotating at flxed ^2 > 0), becomes a stan- 
dard external fleld term in a Heisenberg model. It is more convenient to use 
the Heisenberg models which are natural vectorial reformulations of the Ising 
models for generalized barotropic flow, as shown below. 

The principle of angular momentum conservation and the analogy be- 
tween angular momentum and magnetic moments are both beautifully il- 
lustrated in the famous Einstein-de Haas experiment where a ferromagnetic 
rod is suspended by a thread inside a coil. Upon turning on the current, 
the rod is magnetized, that is, its microscopic magnetic spins are aligned. 
What is surprising is the fact that the rod rotates inside the coil because the 
macroscopic alignment of magnetic moments results in a nonzero net angular 
momentum inside the rod due to large numbers of aligned orbiting electrons, 
and since no torque is applied to the system, the rod reacts by rotating the 
opposite way to conserve angular momentum. 

A related process, albeit one that involves phase transitions, is shown 
in this paper by obtaining exact solutions for the partition functions of the 
spherical Heisenberg models which are the above Heisenberg models plus 
the spherical constraint from flxing the relative enstrophy of the flow. At 
sufficiently hot negative statistical temperatures (with small absolute val- 
ues), we show that the spherical Heisenberg model for generalized barotropic 
fiows goes through a second order phase transition between disordered states 
of local spins (vorticity) at low energy and a global ordered state at very 
high energy where the local spins sum to a total magnetization (net angular 
momentum) aligned with the rotation axis of the solid sphere. 

Reverting to the actual situation where the solid sphere has finite mass, 
this phase transition means that the two component fiuid -sphere system 
undergoes a process very similar to the Einstein-de Haas phenomenon: in 
the high energy ordered state the barotropic fluid layer acquires positive net 
angular momentum (super-rotates relative to the solid sphere) due to macro- 
scopic alignment of local vorticity which compares with the orbiting electrons 
in the ferromagnetic rod acquiring net angular momentum from alignment 
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of magnetic moments; the solid sphere slows its rotation rate to conserve 
total angular momentum, just as the magnetized rod rotates the opposite 
way. Sattinger summarized the phenomenon discussed here succintly in the 
phrase "...many little spins turn into a big spin .." [13] 



4 Heisenberg Model for Barotropic Statistics 

Recall that in the spherical Ising model for barotropic flow, given fixed 
mesh points Xk on S"^ and the voronoi cells based on this mesh [?], we approx- 
imate the relative vorticity by discretizing the vorticity field as a piecewise 
constant function, 

N 

i=i 

where sj = uj{xj) and Hj{x) is the characteristic function for the domain Dj, 
that is 



Hj{x) 



1 X e Dj 
otherwise. 



There is however a more natural vectorial formulation that leads to a 
Heisenberg model for barotropic flows on a massive sphere. Instead of rep- 
resenting the local relative vorticity uj{xj) at lattice site xj by a scalar sj, it 
is natural to represent it by the vector 



Sj — SjTlj 



where fij denotes the outward unit normal to the sphere 5*^ at Xj. Similarly, 
we represent the spin > of the rotating frame by the vector 

h = — Qn 
N 

where n is the outward unit normal at the north pole of S*^. Denoting by 
jji^ the angle subtended at the center of S"^ by the lattice sites Xj and Xk, 
we obtain the following Heisenberg model for the total (fixed frame) kinetic 
energy of a barotropic flow in terms of a rotating frame at spin rate Q, 

^ N N 
= -ifZ^JjkSj-Sk + h-^Sj (1) 
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where the interaction matrix is now given by the infinite range 



167r^ ln(l - cos'^jk) 
^ ^ cos 7,, ^ 



the dot denotes the inner product in and h denotes a fixed external field. 

The Kac-Berlin method [TU] can be modified [2| to treat the spherical 
Heisenberg model which consists of and the spherical or relative enstro- 
phy constraint, 

In addition, Stokes theorem implies that it is natural to treat only the case 
of zero circulation. 

Looking ahead, we note the important fact that the following vectorial sum 
or magnetization 

A N 

i=i 

will turn out to be a natural order parameter for the statistics of barotropic 
fiows on a rotating sphere. 

Can the Heisenberg model H'^ on S'^ support phase transitions? More 
precisely we check what the Mermin- Wagner theorem has to say about : 
(1) it has spatial dimension d = 2, (2) it has a continuous symmetry group, 
namely for each element g G 5*0(3), 

and (3) it has infinite range interaction, that is, for any sequence of uniform 
lattices of N sites on S^, 

An 



lim — Jjk = —CO. 

iT^no AT ^ ^ 



Properties (1) and (2) by themselves would have implied via the Mermin- 
Wagner theorem, that does not support phase transitions, since all d < 2, 
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finite range models with a continuous symmetry group cannot have them. 
However, property (3) violates the finite range condition of this theorem. 
Hence, on S*^ can in principle have phase transitions in the thermody- 
namic limit. 

It is interesting to compare this Heisenberg model Hf^ on 5*^ with the 
Ising type model i^Tv for the same barotropic flow in the last section. The 
Mermin- Wagner theorem there allows H]\f to have phase transitions in the 
thermodynamic limit for a different reason: the Ising type interaction Jjk = 
ln(l — cos7j^) has finite range instead of infinite range but does not 
have a continuous symmetry group, only the discrete symmetry Z2. 

Careful monte-Carlo simulations of this model show that there is one 
negative critical temperature in this model, < for all values of the 
spin rate Q |7j. An extension of the Kac-Berlin method to the spherical 
Heisenberg model for Barotropic fiows on a rotating sphere will show that 
BEC transitions through a symmetry-breaking Goldstone mode to the single 
ground state ip^^, occurs for sufficiently high kinetic energies or small negative 
values of the temperature T. 

5 Solution of the spherical Heisenberg model 
for Q > 

The family of Heisenberg models derived above for the barotropic fiuid 
- solid sphere system in a frame rotating at angular velocity f2 > have 
external fields h{N) = and infinite range interactions 

J ^ 167r^ ln(l - cos7jfc) 
cos7^.fc 

Combining it with the vectorial spherical constraint, 

4:71 ^ 

we obtain an extension of Kac's spherical model [TU] to a one-parameter 
family of spherical Heisenberg models which is parametrized by the size 
of the Voronoi lattice on S"^. 
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This family of spherical Heisenberg models for barotropic vortex statistics 
allows us to model the thermal interactions between local relative vorticity 
ijj{x) and a kinetic energy reservoir at any fixed temperature T. The spherical 
constraint enforces the microcanonically fixed relative enstrophy Q > but 
allows angular momentum in each of the three principal directions to change. 
Similar to the equilibrium condensation process found in the case Q = for 
the spherical Ising model [Lim06a], kinetic energy of barotropic fiow settles 
into a Goldstone symmetry-breaking ground state at very small negative 
temperatures Tc < T < (associated with extremely large energies). Unlike 
the = 0, there is no 3-fold degeneracy in the Goldstone modes and only the 
mode ipio which carries angular momentum that is aligned with the rotation 
axis fin, has a large amplitude. 

The exact solution of the spherical Heisenberg models if^ proceeds along 
similar lines to the Kac-Berlin method for the spherical Ising model. In the 
thermodynamic or continuum limit as A'^ oo, the partition function is 
calculated using Laplace's integral form, 




where moreover, the microstate 

s = {si, ...,sn} e R 
satisfies the zero circulation condition 

N 

■ nj = 0. 
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Thus, 



J a—ioo 

^ I CI ^ """^ ("^ " ^ S ^' 

where 

To evaluate the Gaussian integrals in Z^, we expand the relative vorticity 
vectorfield again, this time, in terms of the spherical harmonics, 

l=l,m=—l 

where n{x) is the outward unit normal to at x. We stress that this expan- 
sion need not include ■0oo(^) — ^ because of the zero circulation condition on 

microstates s. 

Solution of the Gaussian integrals requires diagonalizing the interaction in 
in terms of the spherical harmonics {V'imli^i which are natural Fourier 
modes for Laplacian eigenvalue problems on 5"^ with zero circulation: 

oo I 

1=1 m=—l 

N ^ oo I 

j^k 1=1 m=—l 



where the eigenvalues of the Green's function for the Laplace-Beltrami oper- 
ator on 5"^ are 
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and aim are the corresponding amplitudes. Thus, 

i^fc 1=1 m=-l ^ ^ ^ 



and 



ra+ioo 

Yl daim I A>2(tt) / —exp{N 



m=-l 



m=— 1 



dain, I —exp{N 

a-ioo 



\ m=-l 
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f A>2(«)exp(-ivf; ^ [ 



A 

2N' 



, ^7 \ 2 



Q 



where the order of integration of the term exp {^—N Y11^2 Y^m=-i Q^'imj 
be interchanged by choosing Re(77) = a > large enough. 



5.1 Restricted partition function and non-ergodic modes 

Next we write the problem in terms of the restricted partition function 
(q;io,q;i,±i;/3,(5,Q), that is, 



m=-l 

_[_[ Claim / — exp <( iV 



m=— 1 



27ri 



/3 ^ \ ^2 

4iv + g) ^ 



m=— 1 
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// 00 I /■ 
A>2(a)exp -7V^ ^ ( 
\ 1=2 m=-l ^ 
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Because of non-ergodicity of the condensed modes, we should not inte- 
grate over the ordered modes in this problem, namely aim, which are the 
amplitudes of the 3-fold degenerate ground modes V'lm that carry global 
angular momentum. This often used physical argument in the condensed 
matter literature will for the first time be turned into a rigorous proof here. 
A pertinent and important question arises at this point [5]: how many and 
what are the condensed modes in any given spherical model? We will show 
later that only one single class of modes can have nonzero amplitudes in the 
condensed phase of this problem, namely those belonging to the meridional 
wave number / = 1. 

The statistics of the problem are therefore completely determined by the 
restricted partition function Z^(aio, cti,±i'-j f^j Qj Amplitudes aio, ai,±i of 
the ordered modes appear as parameters in this restricted partition function, 
and will have to be evaluated separately. 

Standard Gaussian integration is used to evaluate the last integral, which 
yields, after scaling P'N = /3, 



Jl>^ \ 1=2 m=-l V 



+ 



a 



Im 



nn 

1=2 m=-l 



IT 



Nr] 



+ 



I3'N 
2 



A 



Im 



provided the physically significant Gaussian conditions hold: for / > 2, 



2 Q 



2/(/ + r 

Then the partition function takes the form 



(2) 



1/2 



Zjj{aio,ai^±i;f3,Q,il) oc 



1 

27ri 



a+ioo 



dr] exp < 



2N 1^1=2 l^m \Q 2 "^'"^ ) 



where the free energy per site evaluated at the most probable macrostate is 
-j.F(?7(/3'),Q,/3') with 



F(r/(/?'),Q,/5')=^(/?') 



m=— 1 



— a?^ — —VtCa 
4 im 2 



m=—l 



2N 



1=2 m 



Nr] (3'N 



X 



Im 
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5.2 Planck's theorem, Saddle points and the Themo- 
dynamic limit 

Provided that the saddle point ?7(/3') can be determined at given inverse 
temperature Planck's theorem states that the thermodynamically sta- 
ble (most probable) macrostate is given by the maximum of the expression 
F{r]{l3'),Q, P'). At positive temperatures, the structure of this expression 
where it concerns the ground modes aim, namely. 



x(aio, ai,±i; Q, ^) = Vi(3') 



1 



m=—l 



Im 



a 



Im 



m=— 1 



and the fact that the saddle point ri{P') must be positive, suggests that 
for any positive value of the saddle point, the expression x ci-nd therefore 
F{ri{i3'), Q, (3') is maximized by ^^=-1 "^im ~ f*^^ > when planetary 
spin VL is small, and by aio < for large /?' > when planetary spin Vt is 
large. At negative temperatures, we expect to find a finite critical point 
where the two opposing parts of x ^'^^ balanced. In order to prove that these 
heuristic expectations are valid, we will solve the restricted partition function 
in closed form by the method of steepest descent. 

The saddle point condition gives one equation for the determination of 
four variables r/, aim in terms of inverse temperature (3' and relative enstrophy 



= — 



OF 
dr] 



-± 



a 



Im 



2NQ 



EE 



1=2 



r0) 

Q 



-r —Mm 



(3) 



where r] = rj{[3') is taken to be the value of the saddle point. Note that 
it does not depend on the planetary spin rate f2 > 0. We note in passing 
that the same equation holds in the = case. There are two natural 
subcases for the saddle point condition, namely, (A) the disordered phase 
(for \T'\ ^ 1) where equation ([3]) has finite solution r]{[3') > 0, and aim = 
for m = —1,0,1; and (B) the ordered or condensed phase (for \T'\ <^ 1) 
where equation ([3]) has finite solution r/(/9') > only when aim 7^ for 
some m. In case (A) which will be solved below, there is no need to invoke 
additional equations of state as the amplitudes aim = for m = — 1, 0, 1. 

Case (B) requires three more conditions to determine the three amplitudes 
aim and the saddle point ?7(/3') > 0. They are provided by equations of state 
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(or Planck's theorem) for the condensed phase (which do not hold in the 
disordered phase): 

dF f2r0) 3'\ (3' 
— h — am -! 



dF (2ri{B') fl'\ 

Thus, a coupled system of four algebraic equations (Ej), (jlj), ([5]) determines 
four unknowns in terms of the planetary spin f2 > 0, the relative enstrophy 
Q > and the scaled inverse temperature The last two equations of state 
for ai^±i implies that either 

ai,±i = Oor [-^ + ^)=0. 

The first equation of state differs from the other two; this represents reduction 
of the S0{3) symmetry that existed in the f2 = case to symmetry in 
the case of nonzero planetary spin. Together these three equations of state 
imply that when > 0, the only possible solution is without tilt, 



(6) 



ai,±i = 0. 



These values of aim will be substituted back into the saddle point condition 
([3]) to yield a single equation that will be solved below. 
The Gaussian conditions ([2]) imply that for / > 1, 



2/(/ + l) Q 

The critical temperature can be obtained from the saddle point condition: 
(A) in the disordered phase at large |T|, 

where the large limit on the RHS is well-defined and finite for any finite 
1/3' I provided 

> r/* = ^ > 0, (8) 
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because then, each term (Z > 2) in the sum is majorized: for negative tem- 
peratures, 

m + " < (-1 + " < (-1) " 

V Q 2l{l + l)J -\ 4^21(1 + 1)) -\ QJ ' 

and for positive temperatures, 

V Q 2l{l + l)J -\A 2l{l + l)J -\QJ ' 

and the corresponding expressions have well-defined positive limits, i.e., for 
all negative and finite /?', 

and for all positive and finite j3' , 

And (B) in the ordered phase at small |T|, 

where a similar argument proves that the RHS is well-defined and finite 
provided ri{P') > rj*. 

This proves that the thermodynamic or continuum limit of the spherical 
Heisenberg model is well-defined for all negative temperatures because 
it turns out (and is shown below) that the saddle point satisfies (|H]) for 
the disordered as well as the ordered phases. Later we will show that this 
thermodynamic limit exists for all positive temperatures as well. 

The large |T| or small \f3'\ saddle point condition in case (A), 

can be solved and has the property that ri{P') \ 1 as — > 0. 

In case (B), when is large, we discuss (i) (3' < and (ii) P' > 
separately. 
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5.3 Negative critical temperature 

For case (i) /5' < 0, a point is reached at (3'^ where 

such that for (3' < (3'^{Q) < 0, 



hm —j^ — y^y^|-- + A;m| < 1. 



iv-00 (3'NQ 

(We note the significant fact that T^{Q) depends hnearly on the relative 
enstrophy Q but does not depend on Q.) In other words, the extreme saddle 
point 

P'Q 

is no longer adequate to solve (JTj) for P' < (3'^{Q) < 0; any larger value rj > rj* 
does not work either because then 

for P' < P'c{Q) < 0. It remains to check that rj < rj* cannot be used. This is 
due to the fact that 

which implies 

1 ^ ' / P' 

iirj<rj* and (5' < (3'^ < 0. 

From this discussion of (i) < < 0, and after substituting the nonzero 
solution ([6]) of the equations of state back into the saddle point equation, 



1=2 m 
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we derive a single equation, 



IQQ V Q ) ^™2iVQ^^V Q ^2/(/ + l)^ 

(13) 

for the saddle point r/(/?', Q, Q) > r/* when < f3'^ < 0. 

The RHS of this equation can be made larger (resp. smaller) than 1 by 
choosing ri{P') < rj* (> rj* resp.) and since Q = J2i=i J2m ^fm the relative 
enstrophy, we must have 

< (l - ia.f„) < 1 

which means that its LHS lie between and 1. Thus, by choosing a suitable 
viP') ^ V* we should be able to satisfy (fT3l) for P' < P'^ < 0. It remains 
to check that this is consistent with the property afg < Q of the ordered 
solution that is, for all fl, 

f v{P',n,Q) P'Y" ^ , .... 

IQQ \ Q ^ a) - - ^ ' 

Thus, to prove that the saddle point condition f|T3l) has solutions ri{j3', Q, Q) > 
rj* for all n > 0, Q > 0, and for all (3' < p'^{Q) < 0, it is sufficient 
to note that for fixed /3' < P'dQ), its RHS(?7(/9')) < 1, decreases as rj 
> rj* increases, while its LHS(r7(/5')) < 1, increases as ?7 > r^* increases; 
and in such a way that RliS{ri{j3')) is surjective on the interval (0, 1) with 
hm^(^f3>)^^RHS{ri{p')) = for any fixed /?' < and RHS{ri*{p'J) = 1, 
and LHS(?7(/3')) is surjective on (0,1) with \im^(^p'-^^^ LHS{ri{/3')) = 1 for 
any fixed P' < /3^,and LHS{f]{(3')) = for the solution 

^yP) = -A :r>v 



of 

IQQ \ Q '4 

Condition (HM then implies that for all > 0, Q > 0, and for all (3' < 
P'ciQ) < 0,the saddle point ?7(/3', Q, Q) satisfies 

r]{(3,n,Q) > — > r] 
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which proves the reflection property in the following remark. 
Remark 1: Since the extreme saddle point 

ri = 

' 4 

satisfies the saddle point conditions ^ and ( flOj) only at the single value of the 
temperature < that separates the disordered phase from the condensed 
phase, but not at other T < 0, we have shown that the usual phenomenon 
known as, sticking of the saddle point in the ordered phase, does not hold 
here. A more appropriate label for this new saddle point behaviour seen in 
the spherical Heisenberg models for barotropic flows on a rotating sphere, 
is jumping and reflection of the saddle point at the negative critical point. 
Indeed the proof above shows that, for all Q > and Q > 0, and for all 
/3' < /3'^iQ) < 0, the saddle pomt t]{/3') > -SlI}cVQ _ ^ > 
We summarize the above results in the physical theorem: 
Theorem 1: (A) For all fl > and Q > 0, the quantity 

^ 1=2 m=-l ^ ' 

has a well-defined and finite limit, called the critical inverse temperature, 

/3',(Q)= hm P'^{Q,N)>-oo, 

that is independent of the rate of spin Q. 

(B) Moreover, the thermodynamic limit exists for the spherical Heisenberg 
models in the sense that for any Q > and Q > 0, the saddle point 
conditions, 



J™ 2iVg^5(, Q ^21(1 + 1] 



are well-defined and finite, and the saddle point satisfies the condition 
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for all (3' <0. 

(C) For all VL> Q and Q > 0, and for all /3' < /3'^{Q) < 0, the ordered 
phase takes the form of the tiltless ( ai^±i = 0) ground mode a;io(/5', Q, Q)ipio 
with amplitude 

which implies that it is aligned with the rotation Q > (super-rotating) and 
is linear in Q. 



5.4 Positive temperature 

For case (ii) /5' > 0, we note that the Gaussian conditions ([2]) are automat- 
ically satisfied since ri{P',Q)/Q > is required of the saddle point of the 
equation 

We deduce from the saddle point condition (ITSl) . that for any Q > 0, and 
any positive l3'{Q) < oo, there is a saddle point rj{(3' ,Q) > associated 
with the disordered phase aio = 0. Otherwise, there is a finite critical point 
P'cciQ) > th^t satisfies the equation, 

I Vn I 

^ 1=2 m=-l 

which is a contradiction since the sum on the RHS does not converge. The 
extreme case, namely r]{f3') = 0, for the saddle point condition holds at 
precisely one point, that is, for /?' = oo. 

We will show next that fllSp has more than one saddle points at all positive 
temperatures. In addition to the disordered phase solution found above, the 
pure ground mode phase is a saddle point ri',aiQ < of f|T5|) . Using the 
solution ([6]) of the equation of state for amplitude aio in f|T5|) gives us the 
final form of the saddle point condition, 

\ leg V Q "^J J ^™2ivg^4rV Q m + i)) 

(16) 
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The alternate saddle point - if it exists - satisfies 



v'i/3',n,Q)>r^if3',Q) (17) 
because the LHS f|T6|) must satisfy the pair of inequalities, 

A useful condition that is equivalent to the upper bound is 



From this we deduce that when the planetary spin 
the LHS( fT6l) can be made to satisfy the lower bound in ( fTSl) by choosing r]' > 
ri'^{P',n,Q) > where LHSdlEl) equals zero at ri'^{P',n,Q) < 00. LHSdlE]) 
equals one at rj' = 00. U Q = Qc, then clearly, ri'^{P', fi, Q) = 0. But for Q < 
flc, any rj' > will satisfy the bounds in (fTSl) . 

Q > flc '■ RHS (fT6|) equals one at the disordered phase saddle point ri{P' , Q) 
which is independent of Q but equals R{r]'^{P' ,fl,Q)) at ri'^{P' ,Q,Q) > 0. 
The RHS( fT6il decreases to zero from the value R{r]'^{P', Q, Q)) while LHS (|T6il 
increases to one from zero as rj' > ri'^{/3' ,Q,Q) increases towards 00. There is 
an alternate (pure ground mode) saddle point solution of f|T6|) since LHSffT6l) 
equals RHS f|T6l) for some rj' G {rj'^^P' ,Q,Q), 00) by continuity. 

On the other hand, if the planetary spin Q is smaller than the critical 
value f2c = \/QJC^, then at 

v'if3\n,Q)=r]{(3',Q), 

the RHS([I6D equals one and the LHS([T6]) equals L^{ri{p' ,Q),n) E (0,1). 
RHS( fT6l) decreases from one to zero as rj' increases from rj{(3' ,Q) to 00. 
LHS (|T6l) increases from Lc{ri{l3' ,Q),Q) towards one as rj' increases from 
ri{P' , Q) to 00. There is again an alternate (pure ground mode) saddle point 
solution of ^ since LHS([I6]) equals RHS([I6]) for some rj' G {ri{(3\ Q) , Q) , 
00) by continuity. 

This also shows that the thermodynamic limit exists for the spherical 
Heisenberg models for all positive temperatures, in the sense that for 
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any Q > 0, f2 > and all (3' > 0, the saddle point condition (IT^ is well- 
defined and finite along the saddle points r]{j3', Q) and rj\l3\ f2, Q). 

It remains to show that the disordered phase is preferred at high positive 
temperatures and the pure ground mode phase with counter-rotation aio < 
is preferred at low positive temperatures. To compare the free energy per site 
of these two phases we ignore for the moment the infinite sum of logarithmic 
terms in F and focus on the part which depends on aio : 



x(aio;/5,Q,^^) = Tl' 



^^^^1' Q + T 



4 



16 V Q ^4 



13' ^fC^ fr]'iP' ,n,Q) 

Q 



+ 



16 V Q 4j \Q 4j 8 V Q 4^ 

16 V Q "^J ' 

The same quantity for the disordered phase is given by 

Comparing them we get the following inequality which implies that the pure 
ground mode phase is preferred: 

16Q ^ Q [ Q + 4 J • ^^"^ 

Fixing > 0, we deduce from that ([20]) holds only if (3' > 

/3',,(fi, Q) where 

16Q Q \ Q 4. J 
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where such a positive value P'^^ < oo exists by virtue of the mean value 
theorem because for /?' near zero, 

iPr^'C v{p',Q) ( v'{p',n,Q) n 

16Q Q \ Q 4 J 

since ri{P',Q) increases as P' decreases, and for P' very large, 

leg Q \ Q ^) 

since ry(/9', Q) decreases down to zero as /3' increases to oo. We used the fact 
that both saddle points Q) and rj\P' , f2, Q) are smooth functions of /3' 
in the range (0, oo). 

Returning to the infinite sum of logarithmic terms in F, 

1=2 m=-l ^ ^ ^ 

we note that the value of this convergent sum for saddle point Q) is big- 
ger than that for ?7'(/3', fi, Q) in view of (|T71) but this difference is logarithmic 
in the difference r]'{P' ,Q,Q) — rj{P\Q) > 0, and is therefore dominated by 
the algebraic difference 

x(aio; Q, ^) - x("io = 0; /5, Q) 

discussed above. 

This completes the proof that the disordered phase is preferred at high 
positive temperatures but the ordered phase is preferred at low enough 
temperatures where, unlike the negative critical point, the threshold value 
P'cciP'i Q) depends on both relative enstrophy Q and planetary spin VL. From 
( IT9|l we deduce that r]'{P' ,Q,Q) is linear in Q. Since i]{P',Q) does not de- 
pend on Q, this implies P'cd^, Q) decreases as planetary spin VL increases. 
Thus, as VL decreases to zero, the threshold value /?cc(^) Q) tends to oo, and 
the disordered phase is preferred at all positive temperatures in the case of 
a non-rotating massive sphere. 

Unlike the critical phenomenology of the barotropic fluid - sphere system 
at very high energies (negative temperatures) which we have shown arises 
from the reflection of the saddle point at the extreme value rj* (the disor- 
dered phase does not satisfy the saddle point condition at negative T' when 
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\T'\ -C 1), its critical phenomenology at positive temperature is not so much 
based on the breakdown of the saddle points as on the system's preference 
for a smaller free energy. Transitions between these positive temperature 
phases for ^2 > are characterized by a greater degree of smoothness than 
its negative temperature counterpart since the free energy is automatically 
continuous at Q) > 0. 

6 Proof of gapless and unique ground mode 
condensation 

An important result needed in the above exact solution of the spherical 
Heisenberg model is the number and type of modes in the condensed phase. 
We will prove that there is only one thermodynamically stable class of modes 
(Z = 1) that has nonzero energy in the condensed phase and thence, exactly 
these three modes are non-ergodic in this problem. The proof is based on 
the existence of multiple saddle points r]'{(5') ^ vW') application of 

Planck's theorem. 

The first two parts of the proof are relatively short. The complete proof 
will be given in the order: (1) at positive temperatures, the only nontrivial 
condensed mode is associated with ccio < signifying a counter rotating solid 
body flow, (2) at negative temperatures, there cannot be gapped nontrivial 
condensed modes where gapped means that there are ergodic modes / in 
between condensed modes /', and (3) at negative temperatures, the only 
relevant nontrivial condensed mode is the ground mode aio > which is 
associated with super-rotating sohd-body flows. Part (3) is longer because 
of the important property that there is more than one saddle points at some 
negative T, namely (i) the pure ground mode saddle point ctio > and (ii) 
the condensed phase «2m 7^ 0, aio < 0. We will show that saddle point 
(i) has higher free energy per site than saddle point (ii) for all values of 
T < Tc < where the latter is condensed. By the extension of Planck's 
theorem to negative temperatures, the preferred macrostate is the one with 
highest free energy, namely the pure ground mode. 
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6.1 (1) condensed modes at positive T cannot have 
wavenumber / > 1 

Assuming that there is a set of condensed modes with single Z > 1 at positive 
T, the finite part of the per site free energy expression F - after dropping the 
infinite sum - has the form 



1 - 



Q 



At 



m=—l m=—l 



E 

m=—l 



a 



Im' 



Similar to the approach in the previous section, the amplitudes aim of these 
modes are fixed by additional equations of state 



dF 



dx 



+ 



Q l{l + 1) 



Oilr. 



0. 



daim daim 
In order for at least one aim to be nonzero, 

Q + 

which contradicts the positivity of /3' and The distinguished ground 

mode CKio can be nonzero because of the inhomogeneous term — yJlC in the 
corresponding equation of state for the amplitude. 



6.2 (2) condensed modes must be gapless at negative 

T 

We will prove next that at negative temperatures, there cannot be any 
nonzero condensed modes /' > 2 that is separated by ergodic modes in be- 
tween. The situation is clear from the simplest case of a single additional 
nonzero condensed mode with I' > 2. Then, / = 2 corresponds to ergodic 
modes which have to be integrated in the Gaussian integrals, but similar to 
the above analysis of the equations of state, a nonzero amplitude ai'm 
implies that + ^ = which in turn means that + ^ < 0, 

violating the Gaussian integrability condition. 
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6.3 (3) the only condensed mode at negative T is aio > 



Assume (for reductio ad absurdum) that more than one class of modes have 
nonzero amphtude in the condensed phase of this problem, namely those 
belonging to wavcnumbers li^h and also / = 1 (this being the distinguished 
ground mode of the problem, has to be part of the condensed phase) where 
^2 = ^1 + 1 = 3. Other cases where I21 h and / = 1 are not consecutive, 
cannot appear in the condensed phase of a solvable spherical model, because 
the Gaussian solvability conditions 



+ ^>o 



Q 2/(/ + l) Q 

are violated as shown above. 

Then, the restricted partition function is given by 

/I Zi I2 

^Im-i ^l\m-i (^hrm 
m=— 1 m=—h m=—l2 



1 ^1 ^2 i>a+ioo 

m=-l m=-h m=-l2 "^""^^ 



N 



a? 

-1 "Zim 



2N ^ 



Q J Z^m=-1 
\ m ^ Q Z^m=-l ^ 



a 



Im 



10 



\ l^l,h,hm=-l ^ V/ y 

Standard Gaussian integration is used to evaluate the last integral, which 
yields, after scaling P'N = j3, 



L 



D{a) exp 



E E 

l^l,ll,l2 m=—l 



2 Q 



TT 



a 



Vn I 

Im \ — W n I JVr; I ^ x 

l=/=l,h,l2m=-l \ Q 2 "^l^ 



provided the physically significant Gaussian conditions hold: for / 7^ 1, /i, /2, 



Q 21(1 + 1) Q 



1/2 
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Then the restricted partition function takes the form 



(aim, a^im, OLi^m'-, P, Q, ^) OC — 

Zm 



a+ioo 



df) exp < 



N 



(f3'\m , ri\ y^l 2 
I 2 Q J Z^m=-1 "«im 

^ \^ 2 Qy Z^m=-l "i2m 



where the free energy per site evaluated at the most probable macrostate is 
-1, F{r]{(3'),Q,(3') with 

h 1 ' ' 



\m=—li m=—l2 m=— 1 



a 



Irn 



4 ^ 



Im 



m=— 1 



(3'X. 



Um 



E 



a 



2N 



A 



The saddle point condition gives one equation for the determination of 
the variables r], aim, (^hmi C(hm in terms of inverse temperature P', relative 
enstrophy Q, and the fixed rate of spin Q > of the planetary frame, 



Q 



a 



Im 



1 V- s-(v{P')^P\ 



,m=—U 



m=—l2 



m=— 1 



2NQ ^ ^ \ Q 



where = r;(/3') is taken to be the value of the saddle point. Equations to 
close the system are provided by equations of state (or Planck's theorem) for 
the condensed phase: for Z = 1, 



= 
= 



dF 
daiQ 
dF 



I h — I Ctio 



dai 



±1 



\ Q ' 2 J 2 
f2r,{P') P'\ . 



OC 
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and for I — li, I2, m — —I, 0, I, 



OOil^m \ Q J 

We assume that aim 7^ fo^ some m = — 1, 0, 1. 

The proof by contradiction continues by further assuming that at least 
one of the I' — I11I2 equations of state are satisfied by nonzero amphtudes, 
say ai^m 7^ for some m — —li, 0, li. Then 



which in turn imphes 



and thence, 



C^hm = 0, 

ai,±i = 0, 
aio 7^ 0. 

The first of the / = 1 equations of state then imphes 

a„ = -|!)c(5^ + f)"<0£or^'<0 (21) 

since flC > 0; and vice- versa for (3' > 0. 

The contradiction at negative temperatures is obtained first by showing 
that although the saddle point equation 
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has solution rj\[3'), aio < 0, 7^ in addition to the pure ground mode 
solution r/(/3'), aio > 0, = at large values of the scaled inverse tem- 
perature /5' < < , these counter-rotating solutions have lower values of 
F' = F{r]'{P'),Q, P') than the same expression F for the pure ground mode 
solution. Here f3'^ < is the most negative inverse temperature for which (|22|) 
has unique saddle point, namely the disordered phase, r], aio = 0, 0:2™. = 0. 
This critical point has the same value critical inverse temperature as obtained 
in the last section because the RHS (122!) is the same as its counterpart in 
the pure ground mode case in the thermodynamic limit. For all f3' < f3'^, an 
argument similar to that used in the previous section to prove the existence 
of the pure ground mode saddle point, can be used here to prove the alter- 
native saddle point equation (122]) has two saddle points, namely the sticking 
one, 

v'i/3') = -^f > «io < 0, «2m ^ (23) 
and reflected pure ground mode solution 

viP') >v* = «io > 0, a2m = 0. (24) 

By the extension of Planck's theorem to negative temperatures, solutions 
aio < 0, 7^ at < are not the most probable and statistically stable 
macrostate because they have lower per site free energy. 

Using the sticking property, we deduce that the free energy expression 
Fir]'i/3'),Q,(3') is given by 

F' 
— lim 

which apart from the first terms in the infinite sum that vanish like InN 
for large A^, and the new value ?]'(/?') < ti{P'), has the same form as F{rj{j3'), Q, j3') 
in the pure ground mode case. Equation ( 123|) and aio < imply that the 
first two terms of F' is smaller than those in F and the third term is bigger. 
So for f2 > large enough, since ri'{j3') is independent of fl, we can make 
F' < F. 

We can do more by using the sticking property of r]\(3'). Substituting 
rj'{(3') = — -[^ and aio < from fl?I]) into the first three terms in F', and using 
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g 4 



r]\P') - ^nCa.o 
1 V- fNv' P'N 



a 



10 



2N 



Q 



+ 



lin 



the reflection property of tlie saddle point ri{f3'), namely, y ^g ^ + yj > 0, 
we obtain the required comparison at P' < 0, that is, 

In other words, for all (3' < 0^ the per site free energy expression F{rj{l3')) > 
F'{rj\(5')). Thus, the pure ground mode condensed phase is preferred in the 
thermodynamic limit. 

This completes the proof of the result that for all values of planetary 
spin at negative temperatures, there is exactly one thermodynamically stable 
saddle point, that is the disordered phase, r/, aio = at negative > 
and the pure ground mode condensed phase, ?7, aio > at negative /?' < 

An important consequence of this result is that the alternative condensed 
phase, r]'{P'), aio < 0, a2m 7^ can in principle be thermodynamically stable 
in other geophysical flow problems such as the generalized Shallow Water 
Equations (GSWE) on a massive rotating sphere. 

7 Solution of the spherical Ising model for 

Q = 

We briefly review the exact solution of the spherical Ising model in the special 
case of a flxed frame and refer the reader to the literature for details In 
view of the sections on the role of energy and angular momentum in the 
formulation of a statistical mechanics for the barotropic fluid - solid sphere 
system, one might ask why we need to solve this case when we have already 
solved the general f2 > case in the previous section. This case applies to the 
situation of a non-rotating inflnitely massive solid sphere which behaves like 
three inflnite reservoirs of angular momentum for each of the three principle 
directions. As will be shown below this case differs from the > case 
because without a distinguished axis of rotation, the ground or ordered modes 



32 



in the problem has SO (3) degeneracy (due to the 3-fold degeneracy of the 
/ = 1 spherical harmonics). Symmetry breaking in the phase transition for 
this case is therefore more explicit than in the > case. In fact it is clearly 
an instance of so-called spontaneous symmetry breaking (SSB) in the general 
Ginsburg-Landau formulation of second order phase transitions. Similar to 
the Heisenberg models, the negative critical temperature here is due to the 
anti-ferromagnetic nature of the logarithmic interaction in the energy. 
The partition function for the spherical Ising model has the form 

f f N ^ \ 

Zn(x J D{s) exp {-13Hn{s)) 6 ' j 

where the path-integral is taken over all microstates s with zero circulation. 
In the thermodynamic or continuum limit as iV — > (X), the partition function 
is calculated using Laplace's integral form, 



Zisf oc 



J D{s) exp i-PH^is)) S (q^ - J2 si ■ 

= j D{s) exp {-(3Hn{s)) J dr] exp (^Q^ ~ ^ ' 

Solution of the Gaussian integrals require diagonalizing the interaction in 
ifjv in terms of the spherical harmonics {^;^}^]^ which are natural Fourier 
modes for Laplacian eigenvalue problems on S'^ with zero circulation. Since 
the ordered modes are associated with / = 1. wc do not need to integrate 
over aim and henceforth discuss only the restricted partition function 



Zn{P, Q; Oiio, q;i,±i) oc J Di>2{a) exp i - 



1=1 m=—l 



1 



a+iocs 

drjexp I -qN I 1 




where the eigenvalues of the Green's function for the Laplace-Beltrami oper- 
ator on S'^ are 

Aim = 1(^1^ ly ^ = 1> m = -I, 0, / 
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and aim are the corresponding amplitudes. 

Next we exchange the order of integration, which is allowed provided 
a > is chosen large enough so that the integrand is absolutely convergent, 
and rescale (3'N = (3, 

ZAf(/3,g;aio,ai,±i) oc /" dr^exp IrjN ( 1 - 7T X] «im ) " ^ ^ ^ 

^TTZ Ja-ioo \ \ ^ m=-l J ^ m=-l 



2 



We stress that 

exp ■ 



rj=-l / 



is not part of the integrand since the ground modes are not integrated. More- 
over, by standard procedures for a Gaussian, we explicitly solve the inner 
integral. 



/ I,(a,e.p(-5:t(^ + Oa.L)^riri 



N I f \ 1/2 

TT 



provided the following physically important conditions hold 

^ + r/^ > 0, Z = 2, ViV, m = -Z, 0, I. (25) 

Thus, 

which can be written in steepest descent form. 



Zjv(/?,(5;Q;io,ai,±i) oc / drjexpN 

J a—ioo 



Z (X lim — / driexp{NF{ri,Q,f3')) 

-j^ pa+ioo 

= — dricxp{-P'g{ri,Q,P')) 

^"""^ Ja-ioo 



exp{-P'g{rj{P'),Q,P')) 
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in the thermodynamic Umit as — > oo, where the free energy per site, after 
separating out the 3-fold degenerate ground states '0io,'0i,±i' given by 

F{rj{(3'),Q,(3') with 

\ m=— 1 / m=— 1 



i=2 m ^ ^ ^ 



2N 

and the total free energy 

/ N 
N-*oo y jj 

The saddle point parameter rj = r]{(3') is determined by solving the following 
set of four equations for {rj, aim) in terms of given values of inverse temper- 
ature P' and relative enstrophy Q. 
The saddle point condition is 

A set of three additional conditions to close the system is given by the equa- 
tions of state for m = —1, 0, 1, 

The last 3 equations have solutions 

Oiim = or — — h P'Xim — 0, for each m. 

This means that in order to have nonzero amplitudes in at least one of 
the ground / condensed states (which are the only ones to have angular 
momentum), — which implies that the inverse temperature must be 
negative, /3' < 0. 
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The Gaussian condition fl25p on the modes with 1 = 2, 



12 2 ' 

can only be satisfied by < when there is any energy in the angular 
momentum containing ground modes. 

Substituting this nonzero solution into the saddle point equation yields 




^ 1=2 m=-l 



c 



where the critical inverse temperature is negative, finite, and inversely pro- 
portional to the relative enstrophy Q, 

' 1--' 



/ — 2 ^Y^ — / 

The saddle point equation provides a way to compute the equilibrium am- 
plitudes of the ground modes for temperatures hotter than < 0, that 
is, 

forT:<T'<0, ^'^rniT') = §^(l-pj. 

m=~l ^ ^ 

The above argument shows that at positive temperatures (low barotropic 
energy), there cannot be any energy in the solid-body rotating modes. In 
other words, there is no phase transition at positive temperatures when 
f2 = 0. This is the spin-lattice representation of the self-organization of 
barotropic energy into a large-scale coherent fiow at very high energies in 
the form of symmetry-breaking Goldstone modes. The reader should com- 
pare the predictions of the spherical model for barotropic vortex statistics 
contained in these formulae with the results of Monte-Carlo simulations in 
the next chapter. In particular, the linear dependence of the negative critical 
temperature Tc = Tc{Q) < on the relative enstrophy should be noted in 
the spherical model solution as well as the Monte-Carlo simulations [?]. 



36 



The free energy per site in the thermodynamic hmit (N — > oo) has the 
form 

f(v,Q,P') = F{r){(3'),Q,(3')^ ^u-Ts 

m=-l \ ^ m=-l / 1=2 m ^ ^ 

1=2 m ^ ^ 



IGtt 
and 



—fir]\ Q,T) = lim — y y In — f -- + A,^^ 



i=2 m ^ ^ ^ 

= hm — yy In— f -- + Aim I - lim — yy 1 



N^oo2N^^ 2T' V 2 / iv^oo2iV 

1=2 m ^ 1=2 m 

1 v^v^, N f 1 ' 
-1 + Urn > > In h Xim 

1=2 m 



when Em=-i «?m > and T^<T' < 0. We used 

hm — W 1 = 1 

N^oo 2N 

1=2 m 

as well as the sticking value of the saddle point 

When T' < < 0, the saddle point equation imphes that 

^ I / A -, \ -1 



Um ^ y y (t^^ + -^X,^) ' = 1 (28) 
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giving 



= -rj{T')T' + lim ^ In (iVr,(r')| + ^A,^) 

1=2 m ^ ^ ^ 

because X]m=-i "^im — ^^'i 

^mT'),Q,T) = -;^(r')-ry'(rOr' + lim2^J]J]ln(iVry(r')^ + ^Aj 



= -,(r') - „'(T')T' + lim 5^ E E 1" ('^"l^')!^ + ^^'•») 

+H-.'(nr|^E E(''(nf.^A,„)" 
-'™5^E E-^'".("(n^ + ^A,j 

1=2 m ^ ^ ^ 

= -r7(T') - yi{T')V + Hm ^ J] J] In (A^r7(ro|^ + ^A. 

1=2 m ^ ^ 

+ry'(r')r' - lim ^ a,. (ry(r')§ + ^A..) ' 

1=2 ra ^ ^ ^ 

= -'<(r')+i-^EE>-(^''(nf + ^A,„) 
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In comparing the free energy per site on both sides of T^, we use the fact 
that the saddle point parameter is stuck at the value 77* for all <T' < 0, 
that is, 



levrT' 



to compute 



— + im > > n VXlra] 

167r ^-00 2N 2T'\ 2 

1=2 m ^ ' 



levr JV-.00 2N 



1=2 m 



Im 



f{T' < r^) = -r](nT' + lim^5^5^1n(iV.^(rO^ + ^Aj 

1=2 rn ^ ^ ^ 



fiT- 



+ A;, 



167r N^oo 2N 2T' 

1=2 m c 

where equation (1281) yields the saddle point parameter in terms of tempera- 
ture, t] = TjiT'). We deduce that the free energy per site is the same on both 
sides of the critical point. This means there is no latent heat involved in the 
phase transition at T'^ < and it is therefore, a second order transition. 
Since 



dT-' 



(T, < T') = -l+ lim — Win— f-- + Az™^ 
^ ^ iv^oo 2N 2T' \ 2 I 

1=2 m ^ ^ 



AttT. 



/N I 



QN 



1=2 m=-l 



-1 



K{T^ < T) 



dT 



f 



— lim 

Af— >oo 



(r < T;) = -r^(T0+lim-L5^5^lnA(^^(T')T'| + Aj 

1=2 m ^ ^ ^ 

^ E E (^(^')-g + 



4ArT' 

-V{T') + K{T' < Tl) - L{T' < 



and 



/N I 



-I -1 



^ i=2 m=-Z 



-1 
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the above expressions for the derivative ^/ at both sides of the critical 
temperature Tc become 



in) = 



hm 

N^oo QAtiN 



'N I 



=2 m=-l 



EE f^--^)' +^(^c) 



Q 



1=2 m=-l ^ 



dT-' 



i67rr; 



/N I 



1™ A AT 



E (^'--9 



=2 m=-i 



+ K{T'J - L{Z) 



where 



< K(T') = hm —S^y^ln^ (-- + Xi^^ < 

1=2 m ^ ^ 



00 



—00 < 



=2 m 



A/™, < 



smce 



-00 



< hm — > > A,^ - - = --^ < 

1=2 m=-l ^ ^ 



and for aU N, 

/]v I 

N 



^EEf-^^--^) <^EEf-^ + ^^-^ ^^-<o- 



1=2 m=—l 1=2 m ^ 

The significant point here is that the specific heat has a discontinuity at 
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T : 



A 



HZ) 



+ 1 



I 



1™ . 7, T 



1=2 m=-l 



-1 



lim 

N-^oo 2N 



1=2 m 



+ ^Im Mr, 



SttT; 1 

H 

Q 2tt 



167rT' Tv^oo AttN ^ ^ \ 2 J 2tx n-.oo 2N ^ -1/2 + lllil + 1) 



lim 

1 1 y^(2/ + l)/(/ + l) 1 1 , ^, 

= hm — > — '-i > hm — > (2/ + 1) > -oo 

27r N^oo N ^ / / + 1 - 2 2ti n~,oo N ' 

1=2 V ' / 1=2 

which is independent of enstrophy Q since L(T^) is and is proportional to 
Q, and which we expect to be positive. 
From the expression 

f(T'^<T') = ^ + K{r)r 

lOTT 

we deduce that /(T^ < T') increases as T' increases away from and also 
that K{T')T' consists of the sum 

-T's(T') + (uiT" ^ 



IGtt 



where s{T') is the entropy per site and u{T') is the internal energy per site. 
At T' = 0, the internal energy per site u{T') = consists entirely of 
energy in the ground modes. At < T' < 0, this represents the fact that 
u{T') — > is that part of the internal energy in the ergodic modes. 

From (125]) . we deduce that rj = t]{T') decreases as T' becomes more 
negative than T^. Then, from 

f{T' < Z) = -,r + lim ^ In (at,^ + ^A.„) 

1=2 m ^ ^ ' 
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we deduce that f{T' < T'^ decreases as T" decreases away from T'^. Here, all 
the internal energy is in the ergodic modes and none in the ground modes, 
and the entropic term —T's{T') becomes more dominant as T' < decreases. 
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